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Abstract. We report on our ^-"-statistic search for white-dwarf binary signals 
in the Mock LISA Data Challenge IB (MLDC1B) . We focus in particular on 
the improvements in our search pipeline since MLDC1, namely refinements in the 
search pipeline and the use of a more accurate detector response (rigid adiabatic 
approximation). The search method employs a hierarchical template-grid based 
exploration of the parameter space, using a coincidence step to distinguish 
between primary ("true") and secondary maxima, followed by a final (multi-TDI) 
"zoom" stage to provide an accurate parameter estimation of the final candidates. 
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1. Introduction 

The Mock LISA Data Challenges (MLDCs) p] have the purpose of encouraging the 
development of LISA data-analysis tools and assessing the technical readiness of the 
community to perform gravitational-wave (GW) astronomy with LISA. The rounds 
so far have been labelled MLDC1 [2], MLDC2 0], and MLDC1B [4]. The challenges 
have consisted of several data-sets containing different types of simulated sources and 
LISA noise. The three types of sources are white-dwarf binaries (WDBs), coalescing 
supermassive black holes (SMBHs) and extreme mass-ratio inspirals (EMRIs). GW 
signals from WDBs will be long-lasting and (quasi-)monochromatic with slowly- 
varyinjU intrinsic frequency /(r); in this sense they belong to the class of continuous 
GWs. In the case of ground-based detectors the typical sources of continuous GWs 
are spinning neutron stars with non-axisymmetric deformations. One of the standard 
tools developed for these searches is the jF-statistic [5], which corresponds to the 
generalized log-likelihood ratio. 

We have applied this method in our MLDC searches, adapting the 
LAL/LALApps [6] search code ComputeFStatistic_v2 used within the LIGO 
Scientific Collaboration to search for periodic GW signals in data from ground-based 



\ In fact, the signals in the MLDCs so far have been strictly monochromatic. Frequency evolution 
is being introduced for the first time in MLDC3. 



Improved T -statistic search for white dwarf binaries in MLDC1B 2 

detectors such as LIGO and GEO 600. We have previously conducted searches for 
WDBs on data from MLDC1 [7J and MLDC2 0(3]. MLDC1B is a rerun of MLDC1 
with different source parameters, and gives us a chance to evaluate improvements in 
our pipeline since MLDC1. Among the issues encountered in our original MLDC1 
analysis were inaccurate determination of a subset of signal parameters due to the use 
of the long-wavelength (LW) limit in modelling the LISA response, and the lack of a 
method to distinguish secondary maxima in parameter space from primary peaks of 
true signals. Both of these aspects have been improved in our MLDC1B pipeline. 



2. Continuous gravitational-wave signals 

A system with an oscillating mass quadrupole moment emits GWs described, far from 
the source, by a metric perturbation h. The WDB signals in the MLDCs 1, 2 and IB 
have been restricted to monochromatic signals with constant intrinsic frequencies /, 
so in an inertial reference frame such as the solar-system barycentre (SSB), the phase 
of this signal can be written as <fi(r) = 27r/(r — T re f), where r rc f is a reference time. 
We refer the reader to [7J for a more complete discussion of the formalism, here we 
only introduce the notation and key results used in the following derivation. The GW 
tensor can be expressed as 

h(r) = A»h^T) , (1) 

where we introduce the convention of an implicit sum J~) ^ over repeated indices /i, v. 
The four amplitude parameters {-4 M } are determined by the overall GW amplitude 
ho, the inclination angle i of the orbital plane, the polarization angle ip, and the initial 
phase (/)q. The explicit relations A^ — A^(ho, l, ip, (j>o) can be found in Eq.(4) of [7J. 
The tensor wave components {h^} depend on the frequency / and the propagation 
direction k of the GW (determined by the sky position of the source), namely 

h x {T)6) = eV(fc) cos0(t) , h 2 (T;8) = *e x (k) cos0(r) , 

h 3 (T;9) = s+(k) sin</>(r) , hi{r;0) = e x {k) sin0(r) , 

where we denote the set of Doppler parameters 9 = {/, k}. The polarization basis £*+,x 
associated with the sky position is defined in terms of the right-handed orthonormal 
basis {£, rj, k} with £ lying in the ecliptic plane and rj in the northern hemisphere, 
namely e+ = £ ® £ — rj <&rj and e* x = £®7? + 77®£. 



3. The ^-statistic method 



The ^"-statistic was originally developed in [5] and extended to the multi-detector case 
in [TU] . A generalization to the full TDI framework for LISA was developed in [TT] . In 
this work we present an approach that unifies the method for space- and ground-based 
detector data, allowing for a more direct application of existing LIGO/GEO600 codes 
to LISA data analysis. 

A "detector" / (here a TDI observable) , provides a linear transformation of the 

tensor metric perturbation /i m (t) into a scalar "signal" h^t) as a function of detector 
time, so the detector output can be written as 

x I (t) = n I {t)+A> l h I fi (t;e), (3) 
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where n 1 it) is the instrumental noise in detector /. Following the notation of [TUlll2j . 
we write the different data-streams x 1 (t) as a vector x(t), and we define the standard 
multi-detector (with uncorrelated noise) scalar product as 

/oo 
x«(f) {S a i{f)]- X y' a {f)df . (4) 

Here we have broken up the observation time into intervals labelled by the 
Fourier-transform of the data in the ath time interval, x* denotes complex conjugation, 
and {S a i(f)} is the two-sided noise power spectral density appropriate to the ath 
time interval. We search for a signal {A S ,9 S } by seeking the parameters {A C ,9 C } 
which maximize the log-likelihood ratio 

L(x; A, 9) = ix\h) - i(h|h) = A^x\h^) - ^{h^h^A" , (5) 

with automatic summation over repeated amplitude indices fi, v. Defining 

Xp{ff) = [x\h p ) , and Mp„(0) = (h„\h u ) , (6) 

we see that L is maximized for given 9 by the amplitude estimator A% = M.^"x u> 
where A4^ v is the inverse matrix of M.^. Thus the detection statistic L, maximized 
over the amplitude parameters A, is T ', where 

2T{x;9)= x^M^ x v . (7) 

This defines the (multi-detector) ^-statistic. One can show that the expectation in 
the perfect-match case 9 = 9 S is E[2T{0 S )} =4 + 14,1, where we used the definition 

\A\ 2 = A^M. liV A v , (8) 

for the norm of a 4- vector A 11 , using M.^ v as a metric on the amplitude-parameter 
space. Note that \A & \ is the (optimal) signal-to-noise ratio (SNR) of the true signal 
{AM- 

4. Modelling the LISA response 

The MLDC data were generated by two different programs: Synthetic LISA [T3] 
simulates a detector output consisting of Doppler shifts of the LISA lasers due to 
relative motion of the spacecraft, while LISA Simulator [TJ] simulates the phase 
differences between laser light following different paths between the spacecraft. In both 
cases the underlying variables are combined with appropriate time shifts to form TDI 
observables which cancel the (otherwise dominating) laser frequency noise [T^l EH [TT] . 
One choice of such TDI quantities is the set of three observables {X, Y, Z}, which were 
used to publish the data of the first and second MLDCs. These observables, which can 
be thought of as representing the output of three virtual "detectors" (which we label 
with the index /), are related to the gravitational wave tensor h through the detector 
"response" , which can be modelled at different levels of accuracy. In the following we 
discuss two such approximations for the response, the simple "long-wavelength limit" 
and the more accurate "rigid adiabatic approximation" . 
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4--1- Long-wavelength limit (LWL) response 

In the LWL approximation the reduced wavelength c/(2nf) is assumed to be large 
compared to the distance L between the spacecraft, which corresponds to a light-travel 
time of T = L/c ^ 17s (assuming equal arm-lengths), and so this approximation 
requires / -C 10 mHz. In this approximation the GW contribution to each observable 
can be modelled as 

^ synthLISA --4T 2 ^ WL :^, (9a) 

x Lis Asim sa _ 2r ^X L: ^ 5 (96) 

where : denotes the contraction of both tensor indices, and d x m = (ri2 ( 8)^2— ^3<8>^3)/2 
is the usual LWL response tensor for a GW interferometer with arms rii and n^. The 
analogous expressions for Y and Z are obtained by cyclic permutations of the indices 
1 — ► 2 — ► 3 — ► 1. In the remainder of this section we will give explicit expressions 
associated with the X variable, with the understanding that the formulas related to 
Y and Z can be constructed by analogy. 

It is convenient to describe the "response" of a gravitational wave detector in the 
frequency domain in terms of a response function i?(/), relating the detector output 
to a "strain" more closely connected to the metric perturbation tensor h, so that 

In the long- wavelength limit, d x w d x and 

to => ' LWL 

i? synthLisA (/) « c hLISA (/) = (i^) 2 ( Ua ) 

This formalism is valid in the regime where the finite lengths of data used to 
approximate the idealized Fourier transforms are short enough that the geometry 
and orientation of the detector doesn't change significantly during this time. 

4-2. Rigid adiabatic approximation 

A more accurate approximation to the TDI response is the so-called rigid adiabatic 
approximation [17] , which for a wave propagating along the unit vector k results in 

where (defining £(fc) = 1 — k • n) 

i27r/Tfcfi/3 

Mf, k) - {e"^ T « « sinc[7r/re(-fc)] + e^^-^ sinc[7r/T£(fc)]} (13) 

is a transfer function associated with the arm along n. Note that this is related to the 
7k(/, k) defined in [T7] by an overall phase, and also that Tr(/, k) reduces to unity in 
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Table 1. Definitions of the long- wavelength (LW), partial rigid adiabatic (pR), 
and full rigid adiabatic (RA) formalisms, in terms of the response function R(f) 
(used to calibrate SFTs) and the detector tensor d 1 (/, A;). 



full name 


label 


response 


detector tensor 


long- wavelength 
partial rigid adiabatic 
full rigid adiabatic 


LW 
pR 
RA 


fi LWL (/) 

R(f) 
R{f) 


d 1 

LWL 

d 1 

LWL 

d ! (f,k) 



the LWL / <C 1/(ttT). For the separation of (fT^|) into a response function R(f) and 
a detector tensor d x (f, k), we choose 

d*(f,k) = k) ^ - X-*(/,£) ^} (U) 

R ( f) e i47rfT 

R(f)= T n m ■ (15) 

sine [2ttj1 ) 

^.5. Calibrated SFTs 

The input to the LAL/LALApps search code consists of Fourier-transformed data 
stretches of duration Tsft, referred to as Short Fourier Transforms (SFTs). This is 
a common data format used within the LIGO Scientific Collaboration for continuous- 
wave searches (e.g., see [IE]). The time baseline Tsft has to be chosen sufficiently 
short such that the noise-floor can be approximated as stationary and the rotation 
and acceleration of the LISA detector can be neglected, and we chose Tsft = 7 days. 

We produce "calibrated SFTs" by Fourier-transforming the raw TDI data and 
applying a frequency-domain response function to produce a Fourier transformed 
strain (including noise) of 

x x (f) = R(f)X(f). (16) 
For our MLDC1 analysis [7] and MLDC2 submission [8] we used the long- wavelength 
approximation -R LWL (/) for calibrating SFTs, but for subsequent analyses we have 
produced "rigid adiabatic" SFTs, which use the full form of R(f) defined in (fl~5|) . 

4-4- Modelling of detector response in different analyses 

Our MLDC1B pipeline includes modifications to implement the full form of d ! (f, k). 
However, a logistically simpler intermediate approximation was also used in the initial 
followup to our MLDC2 work. In this "partial rigid adiabatic" formalism, the more 
precise form of R(f) from (fl~5| is used to construct the SFTs, but the further analysis 
proceeds with the simpler form of d 1 . See table [T] for a summary of the three 
different levels of response approximation considered in this analysis. 

5. Signal templates in the rigid adiabatic formalism 

5.1. Amplitude Modulation coefficients 

In the long- wavelength and partial rigid adiabatic approximations, the strain h (t) at 
the detector is taken to be the contraction of the metric perturbation h{t{r)) with 
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a detector tensor d 1 (t~) which is independent of the frequency and sky direction of 
the signal, but which varies slowly with time due to the change of orientation of the 
detector, in this case as LISA orbits the sun: h 1 ^) = d£ WL (t) : h(t(r)). The template 
"basis functions" ([2|) therefore read as 

h[{t;6) = a{^{t,k) cos0(t(M)), h I 2 {t-6)=b[^{t,k) cos0(t(M)), 
h I 3 (t;0) = a I hwJj (t,k)sin4>(r(t;e)), hl(t;9) = ^ WL (t, k) sM (r(t; 0)) , 

where we have defined the usual amplitude-modulation factors (t, k) = d^ (t) : 

*e+(k) and 6j WL (t, k) = d^ WL (i) : *e x (k). In the full rigid adiabatic analysis, however, 
we need to perform this conversion in the frequency domain, because the response 
tensor d 1 (/, fc) depends on the frequency of the incoming waves. In the rigid adiabatic 
approximation, we model the slow time dependence of the orientation of the detector 
by using a detector tensor d a {f,k) appropriate for the time of the ath SFT. This 
yields the templates 

t(/;^fli(/,fc)c^i(/;9), hi !2 (f;6) = b I a (f,k)^ a (f\0) , (lg) 
< 3 (/;0)=4(/,£)sm71(/;0) , hi A (f;9) = bi(f,k)^fM;0). 

Now the amplitude modulation (AM) coefficients are constructed from the response 
tensor according to 

ai(f,k) = di(f, %):?+(%), &£(/,fc) = d£(/,fe):e x (fe). (19) 
Note that the amplitude-modulation functions a and b now depend on frequency as 
well as sky position, contrary to the LWL case, and they are now complex, due to the 
complex response tensor ([T2]) . 

5.2. Amplitude metric 

Using the rigid adiabatic forms of the templates (fT8|) gives the following form for the 
amplitude metric M.^ defined in ([5]): 

(AC E \ 

C B —E 

—E A C 

\ E C B J 

This form was first exhibited in [Tl . In our notation, the non-zero amplitude metric 
elements for a signal with intrinsic frequency /o are 

a i a i ^ 

Note that the off-block-diagonal element E vanishes in the long-wavelength (and 
partial rigid adiabatic) case because the detector tensor d 1 , an d thus the AM 
coefficients, are real. 

For this project, we have enhanced the ^"-statistic search in the LAL and 
LALApps libraries to allow for AM coefficients which depend on frequency and sky 
direction, and for the more complicated form of the amplitude metric (|20p arising from 
the complex AM coefficients. 



{M MV } 



(20) 
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Stagle-IFO: TDI X 
Wide-parameter .T-'-stat 
Keep N loudest candidates 
above threshold 2T > 2J- t h 



Singlc-IFO: TDI Y 
Wide-parameter ^-stat 
Keep N loudest candidates 
above threshold 2T > 2T th 



Singlc-IFO: TDI Z 
Wide-parameter .F-stat 
Keep N loudest candidate; 
above threshold 'IT > 2J"„, 



Find local maxima within 
metric spheres of mismatch 
m < m LM 



Find local maxima within 
metric spheres of mismatch 
m < m LM 



within metric sphere m c i 



ZOOM in coincidences by 
zoomFactor 



Find local maxima within 
metric spheres of mismatch 
m < m LM 



Keep only maxima coincident 
within metric sphere m r 2 



ZOOM-IN follow-up 
by zoomFactor 

and coherent X ffl (Y — Z) 




Figure 1. The pipeline used for MLDC1B. The pipeline settings used in this 
search were: initial template-grid mismatch mo = 0.25, local-maximum sphere 
of mismatch rriLM = 4.0, first coincidence-stage mismatch m c \ = 0.8, followup 
coincidence mismatch m c 2 = 0.35. The zoom-stage used zoomLevels = 2 and 
zoomFactor =10. 



6. Results and evaluation 

As in the first MLDC, the WDB portion of MLDCfB consisted of seven challenges, 
labelled l.l.la-c and 1.1.2-5. For each challenge, "training" and "blind challenge" 
datasets were provided, generated with different randomly chosen sources. We 
performed the analysis for our challenge entry on the LISA Simulator data, but 
discovered subsequently that inconsistent metadata in this dataset led to a 7.5-second 
time offset and associated systematic error in the initial phase. We are therefore 
presenting here the results of an analysis with the same pipeline on Synthetic LISA 
data. 

6.1. MLDC1B pipeline 

A major limitation of our MLDC1 analysis [7] was the lack of a robust method for 
distinguishing secondary maxima in Doppler parameter space from additional true 
signals with lower .F-statistic values. Starting with our MLDC2 analysis [5], we 
implemented a coincidence condition, where consistent signals are required in searches 
performed will all three TDI variables (X, Y, and Z). This is described in detail in 
a forthcoming paper [9]. Figure [T] illustrates the pipeline used for this analysis. The 
first stage consisted of wide-parameter searches over the Doppler space 9 — {/, k} , 
using each of the three TDI variables X, Y and Z independently. The template 
bank used in this first pipeline stage used an isotropic grid on the sky, with angular 
mesh size da — y/2m / (2n f R OI b/c), with the orbital radius i? or b = 1 AU. The 
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frequency spacing used was df = yJYl mo/(iv T), where T — 1 y is the observation 
time. The mismatch mo used in this first stage template bank was mo — 0.25. The 
expressions for these step sizes were computed from the orbital metric gij |12j . The 
mismatches used in the local-maxima finder and coincidence steps were computed 
from this metric, using the definition m = g^A^A^, where AO 1 are the Doppler 
coordinate differences between two candidates. Local maxima are defined as the 
loudest candidate within a metric sphere of mismatch m < djlm- Coincidence between 
X, Y and Z is defined a having at least one local maximum from each TDI variable 
within a coincidence-window of mismatch m < m c i in the first coincidence step, and 
m C 2 in the second (zoomed) coincidence step. This coincidence scheme was found 
to be effective in eliminating candidates related to secondary maxima. Zooming of 
candidates was achieved by running a search covering 4 neighbouring template points 
in each dimension with a template-grid resolution increased by a factor zoomFactor, 
i.e. m' Q = mo/zoomFactor. The final zoom-stage serves only to increase the parameter- 
estimation accuracy of the final candidates, and is using a coherent TDI combination 
of the noise-independent variables X and Y — Z. 

6.2. Evaluation 

In order to evaluate our errors in parameter estimation, we compare our estimates 
to the injected parameters (provided in the 'key', and denoted by a subscript 's'). 
In Doppler space, Af — f c — / s denotes the frequency error, and s k y is the angle 
between recovered and true sky position; they can be combined into £g = \As\ , 
where eg measures the number of "sigmas" error in the Doppler mismatch between 
the recovered parameters and the injected ones. This is based on the Fisher-matrix 
Tij, which is related to the Doppler metric via = |-4 S | 2 , e.g., see |12j . and 
where m = gijAO 1 A6\ with A9 l = Q\ — 0\. This error measure is chosen in such 
a way that in Gaussian noise it should satisfy E[eg] = 1. The Doppler metric so- 
used here is a simplified "phase metric" ; ideally the full JF-statistic metric should be 
used instead [12] , so this should only be considered a rough estimate of the statistical 
significance of the Doppler errors Ad. For challenges with multiple signals, the Doppler 
mismatch was also used to distinguish found signals (where the candidate parameters 
were within m < 1 of the true signal) from false alarms (candidates having no true 
signal within m < 1). 

The amplitude parameter errors are characterized by 6 a = - — cl 2 , which 

measures the error in the length of the 4D amplitude parameter vector, and c/)a, which 
is the angle between the recovered and true amplitude parameter vectors (calculated 
using the amplitude metric .M M „((9 S ), as in [?])■ These two measures of amplitude- 
error are defined so that for statistical errors due to Gaussian noise, the expectation 
is 2?[5.a] = and the standard deviation of both 5a and 4>a is |-4 S | _1 - The amplitude 
parameter errors can be combined into a "sigma" error, which is = |A^4| /2, with 
A A = Ac — As, defined so that E[ej] = 1. (Note that the use of 6 a for the error in the 
magnitude and €a for the overall amplitude-error is a change in notation from [7j.) 

6.3. Isolated binaries (Challenge 1.1.1) 

This challenge consisted of three separate data sets, each containing one WDB signal at 
an unspecified sky position and within a given frequency band: in 1.1.1a at ~ 1 mHz, in 
1.1.1b at ~ 3 mHz, and in 1.1.1c at ~ 10 mHz. In each case we were able to recover the 
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Figure 2. Recovery of amplitude parameters in Challenges IB. 1.1a (left), IB. 1.1b 
(middle), and IB. 1.1c (right), using LW, pR and RA response models (see tablc[TJ. 
Each plot compares the recovered amplitude 4- vector Ac to the injected signal 
4- vector A a , shown in the plane defined by the two vectors. Gaussian fluctuations 
would lead to a separation of the endpoints of the order |A*4| = tA ~ 2. (Note 
that the two-dimensional projection is slightly misleading: since all vectors do not 
lie in the same plane, the recovered amplitude vectors should be compared only 
with the key and not with each other.) 



Table 2. Recovery of Doppler- and amplitude parameters in Challenges 
lB.l.la-c. The results indicate adequate Doppler-parameter recovery regardless 
of the response model used. As seen previously in figure [2] recovery of amplitude 
parameters improves using the full RA response, especially at higher frequencies. 
Note, however, that the errors e_4 in amplitude parameters with the full RA 
response still appear to be a bit larger than the statistical expectation. 



Challenge 


Rcsp 


A/(nHz) 


</>sk y (mrad) 




S A 


4>A 


M 


IB. 1.1a 


RA 


-0.68 


46.06 


0.54 


0.12 


0.11 


1.14 


(/ = l.lmHz; 


pRA 


-0.68 


61.88 


0.70 


0.12 


0.12 


1.18 


\As\~ 1 =0.07) 


LW 


-0.68 


61.88 


0.70 


0.11 


0.19 


1.57 


IB. 1.1b 


RA 


0.95 


7.71 


0.94 


0.02 


0.10 


1.32 


(/ = 2.9mHz; 


pRA 


0.95 


12.30 


1.03 


0.01 


0.11 


1.37 


\As\~ 1 = 0.04) 


LW 


0.95 


12.30 


1.03 


-0.01 


0.61 


7.40 


IB. 1.1c 


RA 


1.84 


7.49 


0.72 


-0.09 


0.31 


2.37 


(/ = 9.9mHz; 


pRA 


1.84 


5.12 


0.67 


-0.23 


0.25 


2.59 


\Aa\- 1 = 0.06) 


LW 


1.84 


5.12 


0.67 


-0.32 


2.34 


11.60 



injected signal, regardless of the model used for the LISA response, but the recovered 
amplitude parameters were considerably better, especially at high frequencies, using 
the full rigid adiabatic response than the long- wavelength limit. The recovery of 
amplitude and phase parameters is illustrated in figure [5] and tabled) 

6.4- Verification binaries (Challenge 1.1.2) 

In Challenge 1.1.2, the sky position and frequency of 25 "verification binaries" 
was given, while the amplitude parameters of the injected signals were unknown. 
We therefore performed a targeted JF-statistic search at each of the specified sets 
of Doppler parameters, and found the maximum- likelihood estimators A c for the 
amplitude parameters. These results are shown in figure [3] and illustrate the 
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IB. 1.2 Amplitude Parameter Errors 
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Figure 3. Amplitude parameter errors in Challenge IB. 1.2. The quantities 
plotted are denned in section [6.21 We see a continuation of the trends observed 
for isolated binaries: the LW response has phase errors (<Pa) which grow with 
frequency; at higher frequencies, both LW and partial RA exhibit a loss of SNR 
(negative 5^), which is much less pronounced with the full RA response. The 
latter nonetheless displays a slight systematic excess in overall error, i.e. e_4 > 1. 



Table 3. Summary of the number of signals found, signals missed, and false 
alarms in Challenges 1A. 1.3-5 and IB. 1.3-5. We distinguish missed signals with 
\As\ 2 > 40, which should in principle be detectable with our current pipeline, 
and those with \As\ 2 < 40, which are likely to be too weak to pass our 
detection threshold. Challenge IB. 1.3 contained no detectable signal, and the 
only candidate returned by our search was one low-significance false alarm with 
IT ~ 30. 







Found 




Missed 




False 












> 40 


|A 


2 < 40 




Chal 


/ range (mHz) 


1A 


IB 


1A 


IB 


1A 


IB 


1A IB 


1.1.3 


2-7 


15 





5 








20 


1 1 


1.1.4 


3.000 - 3.015 


15 


13 


27 


22 


3 


17 


5 1 


1.1.5 


2.9985 - 3.0015 


3 


3 


30 


29 





12 






performance of the various LISA response models as seen already in Challenge 1.1.1. 
6.5. Recovery of multiple signals (Challenges 1.1.3-5) 

Challenges IB. 1.3-5 were supposed to contain multiple binaries, increasingly crowded 
in Doppler space. Unfortunately, IB. 1.3 was generated with no detectable signals; even 
the loudest had |^4 S | ~ 0.6. To check performance for resolvable, detectable binaries, 
we consider the application of our MLDC1B search pipeline both to MLDC1B data, 
and to the original MLDC1 datasets. The latter search we refer to MLDC1A, to 
distinguish it from the MLDC1 search we reported in [7]. The results are summarized 
in table [3J The corresponding Doppler parameter recovery for 1A.1.3 is shown in 
figure 2J Note that our signal recovery in challenges 1.1.4 and 1.1.5 is still limited 
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Figure 4. Doppler parameter recovery in 1A.1.3. The quantities plotted are 
defined in section [6.21 Most signals were recovered with errors compatible with 
statistical expectations (eg ~ 1); the few outliers may be due to source confusion. 



by source confusion, where the primary maximum of a weaker signal can be modified 
substantially by interference with the secondary peak-structure of stronger signals. 

7. Conclusions and outlook 

We implemented an ^"-statistic search for white-dwarf binaries in the Mock LISA 
Data Challenge based on the LAL/LALApps code developed to search for spinning 
neutron stars in LIGO/GEO600 data. MLDC1B has given us the opportunity to 
improve our MLDC1 search, replacing the long-wavelength response with the rigid 
adiabatic formalism more appropriate to LISA data analysis. We see that our 
amplitude parameter recovery is markedly improved by this enhancement, especially 
at frequencies above 5 mHz. 

Another difference from our MLDC1 search is an improved discrimination between 
secondary maxima and true signals, which is accomplished using a coincidence 
condition between searches on different TDI variables. This is most relevant in the 
case of a large number of signals, such as in MLDC2, and a separate paper will address 
our findings in detail [9] . 

More work remains to improve our handling of source confusion and deal with 
multiple interfering sources, as illustrated by the limited number of sources recovered 
in challenges 1.1.4 and 1.1.5. 
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